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Abstract. We present a systematic experimental study of geometric and statistical 
properties of thermal plumes in turbulent Rayleigh-Benard convection using the 
thermochromic-liquid-crystal (TLC) technique. The experiments were performed in 
three water-filled cylindrical convection cells with aspect ratios 2, 1, and 0.5 and over 
the Rayleigh-number range 5 x 10 7 < Ra < 10 11 . TLC thermal images of horizontal 
plane cuts at various depths below the top plate were acquired. Three-dimensional 
images of thermal plumes were then reconstructed from the two-dimensional slices of 
the temperature field. The results show that the often-called sheetlike plumes are 
really one-dimensional structures and may be called rodlike plumes. We find that the 
number densities for both sheetlike/rodlike and mushroomlike plumes have power-law 
dependence on Ra with scaling exponents of ~ 0.3, which is close to that between the 
Nusselt number Nu and Ra. This result suggests that it is the plume number that 
primarily determines the scaling exponent of the Nu-Ra scaling relation. The evolution 
of the aspect ratio of sheetlike/rodlike plumes reveals that as Ra increases the plume 
geometry changes from more-elongated to less-elongated. Our study of the plume 
area fraction (fraction of coverage over the surface of the plate) further reveals that 
the increased plume numbers with Ra mainly comes from increased plume emission, 
rather than fragmentation of plumes. In addition, the area, perimeter, and the shape 
complexity of the two-dimensional horizontal cuts of sheetlike/rodlike plumes were 
studied and all are found to obey log-normal distributions. 
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1. Introduction 

As an important class of turbulent flows, turbulent thermal convection is ubiquitous in 
nature, ranging from those in the planets and stars, in the Earth's mantle and its outer 
core, and in the atmosphere and oceans to the convection in heat transport and mixing 
in engineering applications. Turbulent Rayleigh-Benard (RB) convection, a fluid layer 
sandwiched between two parallel plates and heated from below, has long been used as a 
model system to study natural convections [H EJ |3]. Two of the important issues in the 
study of turbulent RB convection are heat transport and coherent structures such as 
thermal plumes. In the first one, one tries to understand how heat transported upwards 
across the fluid layer, which is characterized by the Nusselt number Nu = J/(x A /H), 
depends on the turbulent intensity, which is characterized by the Rayleigh number 
Ra = agAH^/vn. Here, A is the applied temperature difference across the fluid layer, 
g is the gravitational acceleration, and a, u, x an d k are, respectively, the volume 
expansion coefficient, kinematic viscosity, thermal conductivity, and thermal diffusivity 
of the convecting fluid. In many theories and experiments, it is often assumed that there 
is a simple power law relation between Nu and Ra, i.e. Nu ~ Ra^ . In the second issue 
one wants to understand the statistical and geometrical properties of thermal plumes, 
which is a localized thermal structure and has been shown to play a key role in many 
natural phenomena and engineering applications, such as in mantle convection where 
mantle plumes just below Earth's crust are responsible for the formation of volcanoes 
(see, for example, Ref. [HE]), in nuclear explosions and stellar convection where plumes 
dominate both the dynamics and the energy transport (see, for example, Ref. [6]). 
The two issues are not independent of each other. It has been shown recently that for 
turbulent RB system most heat is carried and transported by thermal plumes El E] • 
Accordingly, it remains a challenge to establish a quantitative relationship between 
thermal plumes and heat transport. Although a theoretical effort [10] has recently 
focused on this connection, few such experimental studies have been made. 

Thermal plumes generated by a localized hot/cold spot have a well organized 
structure, which consists of a mushroom cap with sharp temperature gradient and a 
stem that is relatively diffusive [HI [12] . Plumes with such a structure are often referred 
to as mushroomlike plumes, which are also observed in turbulent RBC when viewed from 
the side [131 [HI HS1 El HS1 HZl HH HH El [20]. By extracting plumes from temperature 
time series measured locally near the cell sidewall, Zhou and Xia [211 22] showed that 
the size of mushroomlike plumes obey log-normal statistics. However, the morphology of 
the plumes is totally different when one observes from above (or below): thermal plume 
is extended in one horizontal direction but concentrated in the orthogonal horizontal 
direction [231 H31 US [211 (SSI EH [271 EH [211 EOl EI]. Plumes with such a structure 
are often assumed to have significant vertical extent and thus are called sheetlike 
plumes. Puthenveettil and Arakeri [27] studied near-wall structures in turbulent natural 
convection driven by concentration differences across a membrane and found that the 
plume spacings show a common log-normal probability density function (PDF). Zhou et 
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al. [29] further revealed that both the area and the heat content of sheetlike plumes are 
log-normal distributed. Shishkina and Wagner [32] investigated quantitatively geometric 
properties of sheetlike plumes using direct numerical simulations. However, systematical 
experimental studies of geometrical structures of sheetlike plumes are still missing. 
Recently, Funfschilling et al. [30] suggested that when viewed from above thermal 
plumes near the two plates should be better referred to as linelike plumes, as the vertical 
extent of these structures do not seem to be established and it appears more likely that 
they are one dimensional excitations in the marginally-stable boundary layers. 

As two different configurations of thermal plumes coexist simultaneously in 
turbulent RB system, it is natural to ask how sheetlike and mushroomlike plumes 
transform from each other. By using thermochromic-liquid-crystal (TLC) technique, 
Zhou et al. [29J showed, in a cylindrical cell with unity aspect ratio (T = D/H with D 
and H as the inner diameter and the height of the convection cell, respectively), that 
hot fluids (plumes) move upwards, imping on the top plate from below, then spread 
horizontally along the top plate and form waves or sheetlike plumes. As they travel 
horizontally along the plate, these sheetlike plumes collide with each other or with the 
sidewall, convolute and form swirls. As these swirls are cooler than the bulk fluid, they 
move downwards, merge and cluster together. By the symmetry of the system, the same 
process is expected for the morphological evolution of thermal plumes occurring near 
the bottom plate. However, in Ref. [29] the morphological evolution was visualized 
only at a single Rayleigh number (Ra = 2.0 x 10 9 ), and it is not clear whether this 
process is universal for much higher Ra or for cells with different aspect ratios. The 
understanding of this evolution process is of great importance, as the characterization of 
coherent structures is essential to the understanding of turbulent flows in many systems. 

In the present paper, we report new experiments of the temperature and velocity 
fields measured at varying depths from the top plate over the Rayleigh number range 
5 x 10 7 < Ra < 10 11 and the Prandtl number range 4.1 < Pr = v/k < 5.3 and in 
three water-filled cylindrical sapphire cells with aspect ratios 0.5, 1, and 2. We present 
quantitative results on the relationship between plume number density and Nu and 
results that relate the evolution of the plume morphology to the heat transport. The 
remainder of the paper is organized as follows. We give a detailed description about the 
experimental setup and data analysis method in section 2. In section 3, we study the 
z- and i?a-dependence of plume number and discuss the relationship between thermal 
plumes and the Nu-Ra scaling. Section 4 presents study on the geometric properties 
of sheetlike plumes, which are mainly based on data from the aspect ratio 1 cell. We 
summarize our findings and conclude in section 5. 
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Figure 1. A schematic drawing of the convection cell used in the experiments. The 
top and bottom plates are made of sapphire so that the flow in a horizontal plane can 
be visualized and captured. A: the cooling chamber cover, B: the top sapphire plate, 
C: the Plexiglas sidewall, D: the bottom sapphire plate, E: the heating chamber cover, 
F: thermistors, G: nozzles for cooling or heating water, H: stainless steel rings as part 
of the heating and cooling chambers, I: nozzle for filling fluid into the cell, J: nozzle 
for letting air out of the cell, K: a square-shaped glass jacket. 



2. Experimental setup and procedures 

2.1. The convection cell and the experimental parameters 

The experiments were carried out in three cylindrical sapphire cells [331 121] • The 
schematic diagram of the cells is shown in figure [TJ To capture and study the horizontal 
temperature and velocity fields from top, two sapphire discs (Almaz Optics, Inc) with 
diameter 19.5 cm and thickness 5.0 mm were chosen as the top and bottom plates for 
their good thermal conductivity (35.1 W/mK at 300 K) compared with other transparent 
materials. Two chambers, constructed by stainless steel rings H and Plexiglas discs A 
and E, are used to heat the bottom plate and cool the top plate. Each chamber is 
connected to a separate refrigerated circulator by four nozzles, water flows into the 
chamber by two nozzles in two opposite direction and leaves from the other two nozzles 
perpendicular to the inlets. The sidewall of the cell is a vertical tube made of Plexiglas 
(C) with inner diameter D = 18.5 cm and wall thickness 8 mm, respectively. The 
separations between the top and bottom plates H are 9.3, 18.5, and 37 cm so that the 
aspect ratios of the cells are respectively T = 2, 1, and 1/2. Draining (/) and filling (J) 
tubes are fitted on the Plexiglas tube at a distance of 1.5 cm from the top and bottom 
plates. A square-shaped jacket (K) made of flat glass plates and filled with water is 
fitted to the outside of the sidewall [20], which greatly reduces the spatial variation in 
the intensity of the white lightsheet caused by the curvature of the cylindrical sidewall. 
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Figure 2. The optical setup for thermochromatic liquid crystal visualization and 
measurements: S, halogen lamp; LI, concave mirror; L2 and L3, condensing lenses; 
L4, diverging cylindrical lens. 

Four rubber O-rings (not shown in the figure) are placed between the two sapphire 
plates and the steel rings and between the rings and the two Plexiglas discs to avoid 
fluid leakage. Four thermistors F (model 44031, Omega Engineering Inc.) are used 
to measure the temperature difference between the two sapphire plates. To keep good 
contact between the plates and the thermistors, the thermistors are wrapped by the heat 
transfer compound (HTClOs from Electrolube Limited). It is found that the measured 
relative temperature difference between the two thermistors in the same plate is less 
than 3% of that across the convection cell for both plates and for all Ra investigated. 
This indicates the uniform distribution of the temperature across the horizontal plates. 

2.2. Liquid crystal measurements 

The visualization technique employing thermochromatic liquid crystal (TLC) particles 
has been widely used and documented in fluid visualization experiments and was used 
in the present work to visualize the temperature and velocity fields in horizontal fluid 
layers of varying depth from the cell's top plate. Two types of TLC microspheres 
(Hallcrest, Ltd.) were used in the experiments. One type (model R29C4W) was used 
for low-A experiments with the mean bulk temperature To = 30 °C and Pr = 5.3, 
and the other type (model R40C5W) was used for high-A measurements with To = 42 
°C and Pr = 4.1. Both these particles have a mean diameter of 50 fim and density of 
1.03 — 1.05 g/cm 3 , and were suspended in the convection fluid in very low concentrations 
(about 0.01% by weight), at which the influence of TLC particles on the fluid can be 
neglected. The peak wavelength of light scattered by these particles changes from red 
to green and then to blue within a temperature window of 4 °C from about 29 to 33 
°C for R29C4W, and of 5 °C from about 40 to 45 °C for R40C5W. In the experiments, 
the mean bulk temperature was set to 30 °C for R29C4W (Pr = 5.3) and to 42 °C for 
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Figure 3. Images of TLC microspheres with the camera exposure time of 0.02 s taken 
at 2 mm from the top plate at (a) Ra = 1.1 x 10 9 , (b) 3.0 x 10 9 , and (c) 1.0 x 10 10 . (d), 
(e) and (f), which correspond to (a), (b) and (c), respectively, are extracted sheetlikc 
plumes with background in the original images removed. 



R40C5W (Pr = 4.1), so that the background fluid appears blue and the red and green 
regions correspond to cold fluid, i.e. cold plumes. Figure [2] shows a schematic diagram 
of the optical setup for the experiments. A halogen photo optic lamp (S) with a power 
of 650 W was used as the light source. One concave mirror LI and two condensing lenses 
L2 and L3 were used to collect the light from 5* and focus it onto the central section of 
the cell. A horizontal sheet of white light, generated by a diverging cylindrical lens L4 
and then projected onto an adjustable slit, passed through the cell parallel to the top 
plate. The thickness of the lightsheet inside the cell is approximately 3 mm. A Nikon 
D1X camera, with a resolution 2000 x 1312 pixels and 24 bit dynamic range, was placed 
on the top of the cell to take photographs of the TLC microspheres. With short camera 
exposure time (0.02 s), the captured photographs give the instantaneous temperature 
field, and with long camera exposure time (0.77s) they will in addition show the the 
trajectories of the particles. 

To study thermal plume properties quantitatively, we counted plume numbers and 
extracted sheetlike plumes. For sheetlike plumes, we took 150 to 300 consecutive images 
(at z = 2 mm and with camera exposure time of 0.02 s) at 30 to 60 s intervals for a given 
Ra, so that two successive images in each sequence are statistically independent. A cold 
sheetlike plume is extracted by first manually drawing a contour around its perimeter 




Figure 4. Images of TLC microspheres with the camera exposure time of 0.77 s taken 
at (a) 2 mm and (b) 2 cm from the top plate at Ra = 1.0 x 10 9 . 

and then using a software to collect all pixels enclosed by the contour. A program is 
then used to calculate the perimeter P and the area A of each extracted plume. To 
ensure that plumes are identified correctly the operator uses knowledge gained from 
viewing movies of plume motions. A total of 3000 to 6000 plumes are identified from 
the 150 to 300 images, which are not large numbers but should be indicative of the 
statistical properties. Figures E(d)-(f) show three examples of extracted plumes with 
background removed, which correspond to images of TLC microsperes in figures [3(a)- 
(c), respectively, for three different values of Ra obtained from the unity-aspect-ratio 
cell. For mushroomlike plumes, we took sequences of images at varying depth z from 
the top plate and with camera exposure time of 0.77 s at 30 s intervals for each Ra. 
Each sequence consists 200 images for a given depth. A cold mushroomlike plume was 
identified as an object with nonzero vorticity and with temperature (color) much lower 
than that of the background fluid [29J (see, e.g., figure E(b)). As mushroomlike plumes 
are partly entangled together, they cannot be separated completely. Nevertheless, we 
can identify them individually as swirls with cold temperature and thus count their 
numbers. The cold sheetlike plumes were also identified and counted from these images, 
as objects with a linelike shape and lower temperature. 

3. Plume number statistics 

3.1. Morphological evolution of thermal plumes 

The process of morphological evolution between sheetlike and mushroomlike plumes has 
been revealed by Zhou et al. [29] in a cylindrical cell. Here, the same process was also 
observed for all values of Ra and V investigated in the experiment. Figure 11(a), taken 
at 2 mm from the top plate, shows how this evolution comes about. One sees that near 
the top plate the motion of TLC microspheres appear to emanate from certain regions 
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Figure 5. The process of morphological evolution between sheetlike and mushroom 
plumes. 



or "sources" with bluish color, suggesting that hot fluids (plumes) are moving upwards, 
impinging on and spreading horizontally along the top plate. Along the particle traces, 
the color turns from blue to green and red, implying that the wave fronts are cooled 
down gradually by the top plate (and the top thermal boundary layer) as they spread. 
As they travel along the plate's surface, sheetlike plumes collide with each other or with 
the sidewall. As different plumes carry momenta in different directions, they merge, 
convolute and form swirls (hence generating vorticity). As these swirls are cooler than 
the bulk fluid, they move downwards, merge and cluster together (see, e.g, figure IHb)). 
We note that, by symmetry of the system, the same morphological evolution of thermal 
plumes should occur near the bottom plate. The physical picture of this evolution 
process is illustrated in figure El 



3.2. Depth- dependent properties 

Because of the collision and convolution of sheetlike plumes and the transformation 
from sheetlike to mushroomlike ones, it is expected that the number of sheetlike plumes 
should decrease from the top plate, while that of mushroomlike ones should increase, 
and this is indeed the case as shown in figure El Figures Eta), (c), and (e) show the mean 
numbers of (cold) sheetlike plumes N^ eet versus the normalized height z/5 t h obtained 
from the cells with T = 2, 1, and 0.5, respectively. (Here, 5 t h is the thickness of thermal 
boundary layer, which was measured in [3H [35].) One sees that Npl* eet decays quickly 
away from the top plate and the behavior is similar for all Ra studied. Beyond the 
depth z ~ 155^, one can hardly identify any sheetlike plumes from the acquired images. 
Note that, in addition to transforming into mushroomlike plumes, the rapid reduction 
of N^ eet with depth could also be attributed to the limited vertical extent of sheetlike 
plumes themselves, i.e. plumes, whose vertical extents are shorter than a certain height 
h, could not be captured by the images that are taken farther than h from the top 
plate. We further note that the ^-dependence of plume numbers N^ eet may be well 
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Figure 6. Mean numbers of (a), (c), and (e) (cold) sheetlike plumes N^ eet and (b), 
(d), and (f) (cold) mushroomlike plumes N^ ush as functions of the distance z for all 
measured Ra. (a) and (b) for V = 2, (c) and (d) for T = 1, and (e) and (f) for T — 0.5. 

described by a decreasing exponential function, i.e. Npf eet = Noe~ brZ ^ 5th with fitting 
parameters Nq and b r for nearly all measured Ra. Hence z e = l/b r may be used as 
a typical height of vertical extent of sheetlike plumes. Figure El^a) shows the scaled 
height z e /Sth as a function of Ra. It is seen that z e /5th increases with Ra for the T = 2 
and 0.5 cells, while keeps nearly constant for the T = 1 cell. We can not judge the 
significance of the different behaviors of z e /8 t h varying with Ra in cells with different V . 
However, z e varies only in the range of a few times the thickness of thermal boundary 
layer {25 t h < z e < 55 t h)- This suggests that the vertical extent of the so-called sheetlike 
plumes does not have enough spatial extension in the vertical direction to form sheets. 
For mushroomlike plumes, the mean numbers of (cold) mushroomlike plumes N^ ush 
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Figure 7. (a) z e /5 t h and (b) zc/5th vs i?a for V = 2 (up-triangles), 1 (circles), and 0.5 
(down-triangles). The dashed line in (b) marks the mean value of z c for Ra > 3 x 10 s . 
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as a function of the scaled depth z/H is shown in figures EJ^b), (d), and (f). With 
increasing 2, N^ ush for all i?a first increases rapidly, then decreases and finally remains 
approximately constant for positions z > 0.2H. When ]Sf^ ush first increases with z, it 
crosses over with N$ eet at some crossover-depth zq- The normalized crossover depths 
zq are shown in figure C^b). It is seen that z c ~ 2.55th when Ra < 3 x 10 8 . For 
Ra > 3 x 10 8 , although the data points seem to be somewhat scattered, they exhibit 
no clear dependence on Ra and have a mean value of ~ 4.845^. This implies that 
the transformation occurs within a region, whose height is only associated with the 
thermal boundary layer thickness. If one takes this value as the typical vertical extent 
of sheetlike plumes, these results again suggest that the vertical extent of what is called 
sheetlike plumes is only the order of several 5th, not much larger than it, and thus not 
large enough to have a sheetlike shape. This will be further discussed in section 14.41 
After N^l ush attains its maximum value around z = z p , it then drops sharply beyond 
the depth z ~ 0.2H, because of the mixing, merging, and clustering of mushroomlike 
plumes. As found by Zhou et al. [29], this region corresponds to the region of full width 
at half maximum of vertical vorticity profile and can be used as a quantitative definition 
and measure of the mixing zone [361 EI] • It is further found that the peak position in the 
profiles of N£f ush , z p , has no obvious i?a-dependence for each T, but decreases slightly 
from 0.08H for F = 2 to 0.03F for T = 0.5. 
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Figure 8. i?a-dependency of sheetlike plume number density Npf leet /(irD 2 /4) (dark- 
green symbols), counted just outside the thermal boundary layer, and mushroomlike 
plume numbers N'^ ush / (tt D 2 / A) (red symbols), counted at the peak of the height- 
profile of N^ ush , for the r = 2 (up-triangles), 1 (circles), and 0.5 (down-triangles) 
cells. For both N^ eet and N™ ush , the V = 2 (0.5) data have been shifted down (up), 
so that they agree with corresponding data for Y = 1 in their overlap regions of Ra. 
Solids lines are power-law fittings to the data shown in the same color. The fitted 
scaling exponents are 0.29 ± 0.03 and 0.32 ± 0.03 for N^ eet and N™ ush , respectively. 

3.3. The Rayleigh number dependency 

We next focus on the Ra- dependency of plume number density. For sheetlike plumes, 
the number was counted just outside the thermal boundary layer and the counted 
number, N^ 6 ^, for each cell can be well described by a power law N^ eet / '(irD 2 / '4) = 
a s (T)Ra^ 3 ^ r \ with two fitting parameters a s and f3 s as functions of V. Here, the fitted 
values of a s and f3 s for all three aspect ratios are listed in table [TJ It is seen that a s (T) 
varies with the cell's aspect ratio T while /3 S (T) is approximately the same for all three 
cells. As the range of Ra is limited for each aspect ratio, we shifted the plume number 
densities for the T = 2 and 0.5 cells to agree with data for the T = 1 cell in their 
respective overlapping range of Ra. We then fitted a single power-law to the combined 
data set from the three cells that span almost three decades of Ra. Figure [8] shows the 
shifted numbers, normalized by the area of the top plate irD 2 /4, as well as the plume 
number from the T = 1 cell as dark-green symbols. It is seen that a single power law 

N$ eet /(TTD 2 /4) = a s Ra Ps with a s = 1.4 and 0„ = 0.29 ± 0.03 (1) 

can be used to well describe the data from all three cells. 

For mushroomlike plumes, one sees from figures Mjo), (d), and (f) that there is a 
peak on the height-profile of mushroomlike plume number N£f ush for all three cells and 
for all measured Ra. We use the value at the peak of the height-profile of mushroomlike 
plume number to examine the i?a-dependence of N"l ush . It is found that a power 
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Figure 9. i?a-dependence of the mean aspect ratio of sheetlike plumes (7pf eet ) vs. 
Ra in a log-log plot. 

law relation N^ ush / (n D 2 / 4) = a m (r)RaP m ( r ^ can also be used to describe the in- 
dependence of N£f ush . Again, the fitted (3 m (T) is nearly T-independent (see tabled]). 
Figure [S] shows the .Ra-dependence of N£f ush / '(it D 2 / '4) (red symbols). Here, the T = 2 
and 0.5 data have also been shifted to agree with the V — 1 in their respective overlap 
regions of Ra. It is seen that N^ ush shows similar trend with Ra as N^ eet . Again, the 
relationship between N"{ ush and Ra can be well represented by power-law fits 

N^ ush /(TxD 2 /A) = a m Ra^ m with a m = 0.41 and (3 m = 0.32 ± 0.03. (2) 

To make a quantitative connection between thermal plumes and the heat flux, we 
use the previous experimental finding from both Eularian and Lagrangian measurements 
that the heat flux are mainly carried by thermal plumes [9j. As mushroomlike plumes 
are evolved morphologically from sheetlike plumes [29], the total heat flux carried by 
either type should be the same, i.e., 

Nu ~ N^ e£t F^ ieet ~ jymushpmush _ Ra /3 (3) 

Here, F^ eet and F^ ush are the mean heat flux carried by individual sheetlike and 
mushroomlike plumes, respectively, and ~ 0.3 ± 0.02 in the Ra and Pr ranges of the 
experiment (37] EE]. Recall that the ifa-dependence of Nfi eet [Eq. ©] and {N^ ush ) peak 
[Eq. (T5])] discussed above, one sees that the scaling exponents of the Nu-Ra relation and 
of the i?a-dependence of plume numbers are the same within experimental uncertainty. 
This indicates that the heat flux transported by individual plumes is nearly independent 
of the turbulent intensity, i.e. Ra, and hence the Nu-Ra scaling relation is determined 
primarily by the number of thermal plumes. 

4. Geometric properties of sheetlike plumes 

4-1. Aspect ratio of sheetlike plumes 

A striking feature that one can observe from figure [3] is that, as the Rayleigh number is 
increased, the morphology or geometry of sheetlike plumes change from more elongated 
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in a log-log plot. It is seen that (7pf ee *) decreases with increasing Ra and the relation can 
be described by a power-law with a scaling exponent —0.23. We note that (7~f eet ) can 



Figure 10. i?a-dependence of (a) the mean normalized area (AA/ttD 2 ) and (b) the 
mean normalized perimeter (P/ D) of sheetlike plumes. 



shape to less elongated and more fragmented. The extracted horizontal cuts of sheetlike 
plumes may be characterized by a typical length I and a typical width w. With increasing 
Ra, the length of sheetlike plumes seems to decrease, while the plume's width seems 
to increase. To describe these geometric properties quantitatively, we define the aspect 
ratio of sheetlike plumes, 7pf eet , as 

P = 2(l + w), 

A = Iw, (4) 
. I > w. 

Figure [9] shows the i?a-dependency of the mean aspect ratio of sheetlike plumes (7^ eet ) 

, . ,..,1;,.,, r> \\",, ,.,,-(,, (i,.„ /. 

also be described by a logarithmic function of Ra and we cannot definitively conclude 
which one, between power-law and logarithm, is a better choice. Figures [TOTa) and 
(b) show Ra- dependence of (a) the mean normalized area (4A/nD 2 ) and (b) the mean 
normalized perimeter (P/D) of sheetlike plumes. Again, these two quantities are found 
to decrease with increasing Ra. All these behaviors are well illustrated in figure [31 and 
suggest that when the flow becomes more turbulent, the increased mixing, collision, and 
convolution can more easily fragment large-sized sheetlike plumes into smaller ones. 

Figure [TT1 shows the plume area fraction f p i as a function of the Rayleigh number in 
a log- log plot. It is seen that f p i increases with increasing Ra and a power-law function 
with a scaling exponent 0.23 can be used to describe well the relationship between f pt 
and Ra. Using the shadowgraph technique, [30J studied i?a-dependence of the plume 
area fraction. Although the scaling range is limited, they found a power-law scaling of 
f p i with Ra for f p i < 0.1 and the obtained scaling exponent is around 2. 

Next we examine the relationship between plume perimeter and area. Figure 
[T2Ta) shows a scatter plot of the normalized perimeter P/D and the normalized size 
(AA/nD 2 ) 1 / 2 of sheetlike plumes, which contains a total of 6071 sheetlike plumes, 
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Figure 12. (a) Scatter plot of the normalized perimeter P/D vs the normalized size 
(AA/ttD 2 ) 1 ^ 2 of sheetlike plumes at Ra = 3.0 x 10 9 . (b) The conditional average 
(P\A)/D on the normalized size (AA/ttD 2 ) 1 / 2 for the same data as (a). Inset is the 
plot of fractal dimensions of sheetlike plumes' boundary d A vs Ra. Solid lines in both 
(a) and (b) are power-law fits. 

extracted from a sequence of 260 images at Ra = 3.0 x 10 9 . Here, the perimeter and 
area are normalized by the diameter and the area of the conducting plate, respectively. 
Although the data look somewhat scatter, one sees that the perimeter increases with 
area and all data points can be fitted by a power-low function, 

(P/D) ~ (yj4A/nD 2 ) dA with d A = 1.50. (5) 

To better explore this feature, the conditional average (P\A)/D on the normalized size 
(AA/nD 2 ) 1 / 2 is shown in figure GJ^b). In the fi gure, a good scaling range can be seen 
and the solid line is a power-law fit (P/D) = 28.1 x (AA/nD 2 ) 1 - 5 ^ - 02 . We find that 
such a scaling behavior exists for all Ra investigated. The inset of figure [T27b) shows the 
i?a-dependence of d A . One sees that these d A have a mean value of 1.50, but decrease 
from 1.53 to 1.46 when Ra increases from 6.7 x 10 8 to 6.3 x 10 9 . These measured values 




Figure 13. (a) Scatter plot of the shape complexity SI2 vs the normalized size 
(AA/ttD 2 ) 1 / 2 of sheetlike plumes at Ra = 3.0 x 10 9 . (b) The conditional average 
(l^l-A) on the normalized size {AA/irD 2 ) 1 /' 1 with a solid line as a power-law fit. Inset: 
da vs Ra. (same data as in figure [9]) 



of d,A are further found to be larger than the value of 1.29 found for isoconcentration 
contours of passive scalars measured in the same system [40J . 

Mathematically, as introduced by Mandelbrot [H] and Lovejoy [42J, is the fractal 
dimension of the boundary of sheetlike plumes and satisfies 1 < oIa < 2. Fractal 
dimensions in turbulence have been widely studied (see, for example, Ref. [13]). In the 
present case, our results seem to suggest that the boundary of sheetlike plumes is fractal 
with the dimension of around 1.50 and the slight decrease of cLa may suggest that the 
shape of sheetlike plumes becomes smoother when the flow becomes more turbulent, as 
a result of the increased mixing. However, the extracted sheetlike plumes, as shown in 
figure El do not look like fractal objects. To understand this and to see if applying the 
machinery of fractal analysis to the geometric properties of plume can shed some light 
on the problem, we study the shape complexity of the plumes. 



4-2. Shape complexity of sheetlike plumes 

The geometric complexity of an object can be characterized by its shape complexity, 
which is a dimensionless ratio between its area and volume. The shape complexity was 
previously used to study the shape of isocontours for passive scalar fields in turbulence 
[39]. We introduce it here to describe the shape of sheetlike plumes. For a two- 
dimensional (2D) closed contour, the shape complexity is defined as follows, 

where P and A are the perimeter of the contour and area inclosed by the contour, 
respectively, and the subscript 2 refers to 2D. As a circle has the minimum perimeter 
1\J nA among all 2D objects with the same area A, the shape complexity satisfies 
1 < Q 2 < 00 and can be used to describe the departure of an object from the shape 
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Figure 14. (a) i?a-dependence of the mean shape complexity (f^) of sheetlike plumes. 
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of a circle. For a fractal object, a larger Q 2 implies a rougher contour (or surface) the 
object. For a non-fractal object, a larger Q 2 signifies that the object is more different 
from a circle (e.g. it is more elongated). 

Figure [131(a) shows a scatter plot of the shape complexity Q 2 and the normalized 
dimension (or size) (AA/ttD 2 ) 1 / 2 of sheetlike plumes, obtained from the same data set 
as those in figure [12j It is seen that the mean trend of the relation between the size 
and shape complexity can be captured by an increasing function, i.e., the plume with 
larger size has a higher probability of possessing a larger Q 2 . To illustrate this mean 
trend more clearly, the conditional average (f^A) on the normalized size (AA/11D 2 ) 1 / 2 
is plotted in figure fTBT b). The solid line in the figure is a power-law fit, i.e., 

Vl 2 ~ (^AA/nD 2 ) dn with d Q = 0.50 ± 0.01 (7) 

and the reasonable scaling range can be found. Here, dn ~ 4a — 1 is an immediate 
consequence of the definition of shape complexity [Eq. ([6])] together with the power- 
law relation [Eq. The inset of figure ITBTb) shows the i?a-dependence of d^. One 
sees that da decreases slightly with increasing Ra, which may be understood from the 
decrease of dA- 

Figure [TW a) shows the Ra- dependence of the mean shape complexity (f2 2 ). It 
is seen that (0 2 ) decreases with increasing Ra. This suggests that, because of the 
increased mixing, collision and convolution of sheetlike plumes, the shape of sheetlike 
plumes become closer to the shape of a circle when the flow becomes more turbulent. 
This also implies the decreases of plume's aspect ratio with increasing Ra. Indeed, with 
the definitions of ^ eet [Eq. (0}] and Q 2 [Eq. ©], one can obtain the relation between 
of f2 2 and 7pf eet , i.e., 

1 . n^— 1 . 




Figure 15. PDFs of (a) normalized area 4A/ttD 2 , (b) normalized perimeter P/D, (c) 
shape complexity f2 2 ; (d) normalized length l/D, (e) normalized width w/S t h, and (f) 
aspect ratio sheetlike plumes, all obtained at Ra = 3.0 x 10 9 . Solid curves are 
fittings of log-normal function to the respective data. 



which implies that SI2 decreases with decreasing 7^ * when 7^ > 1. Furthermore, 
f^2 is proportional to 7pf eet when 7pf eei ^> 1. Accordingly, the decrease of (f^), shown 
in figure Ul a). may be understood as the result of the change in aspect ratio of the 
sheetlike plumes, i.e. the decrease of (7pf" eet ) (see figure E]), as 7pf" ee * is much larger than 
1 for most of situations. To see this more clearly, we study the ratio Cr between (fl 2 ) 
and A^{J {lpi eet ) H — /, l heet) )- Note that Eq. (JS} is valid only for an individual plume, 

while Cr is in an average sense. Figure flUf b) shows Cr as a function of Ra. It is seen 
that Cr ~ 1 for nearly all Ra investigated. This result demonstrates clearly that the 
decrease of VL 2 is caused by the decrease of the plume's aspect ratio, rather than by a 
decreased roughness in the case of a fractal object which the sheetlike plumes are not. 



4-3. Distributions of geometric measures of sheetlike plumes 

In a previous study, Zhou et al. [29] found that the area of sheetlike plumes obeys log- 
normal statistics. Here, we further investigate distributions of other geometric measures 
of sheetlike plumes. Figures fT5T a)-(f ) show the measured PDFs of normalized geometric 
measures, i.e. 4A/nD 2 , P/D, Q 2 , l/D, w/5 t h, and 7pf ,eet , of sheetlike plumes. It is seen 
that all these quantities have a same distribution, i.e. the log-normal distribution. The 
same distributions were also found for all of these quantities at other Ra investigated 
in our experiments. Together with the log-normal distributions found for mushroomlike 
plumes [21J, these findings suggest that the log-normal distribution is universal for 
thermal plumes and log-normal statistics may be used to model them, at least in 
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turbulent RBC. In addition, this log-normal statistics of thermal plumes is different 
from that of passive scalars measured in the same system, which is found to obey a 
log-Poisson statistics |i0] . 

4-4- Three-dimensional structures of sheetlike plumes 

Finally, we study the 3D structures of sheetlike plumes. Results obtained in section 
13.21 have suggested that the vertical extent of what is called sheetlike plumes is only 
of the order of the thermal boundary layer thickness, and hence not large enough to 
form a sheet. To see this more clearly , we used a tomographic reconstruction technique 
to construct the 3D image of thermal plumes from sequences of 2D images acquired 
near the top plate of the cell. To achieve this, the convection cell was placed on a 
translational stage and, as the cell traverses continuously at a speed of ~ lcm/ sec, a 
series of photographs of TLC microspheres were recorded by a Nikon D3 camera (3CCD, 
with a resolution of 2397 x 1591 pixels) operating at 11 frames/s. As the speed of the 
the horizontal motion of the sheetlike plumes are about OAcm/sec, the sequences of the 
horizontal slices of the plumes may be regarded as taken at approximately the same 
time. In post-experiment analysis, 2D horizontal cuts of cold plumes were extracted 
from the images taken in each run and a MATLAB script was used to reconstruct 3D 
thermal plumes from these extracted 2D cuts. Figure [ToT a) shows an example of the 
reconstruction of 3D thermal plumes near the top plate. From the figure, one can see 
how mushroom-like plumes are formed by the convolution (or spiraling) of sheetlike 
plumes. Figure [ToT b) shows this process more clearly, which is an enlarged region of 
the image. Another noteworthy feature is that most plumes have a one-dimensional 
structure, rather than a sheetlike shape. This geometric feature may be illustrated 
more clearly by figure floTc). which shows the reconstruction of an enlarged individual 
3D plume. It is seen that both the height and the width of this plume extend only to a 
few millimeters, while its length extends to a few centimeters. This is the most direct 
evidence thus far that shows thermal plumes near the conducting plates are only one- 
dimensional structures with horizontal length being much larger than their horizontal 
width and vertical extent, which has been suggested previously [30]. Therefore, we may 
hereafter call thermal structures near the conducting plates rod-like plumes. Future 
investigations will be focused on the geometric properties of 3D structures of both rod- 
like and mushroom-like thermal plumes. 

5. Conclusions 

In this paper, we have presented a detailed experimental study of the temperature and 
velocity fields in turbulent Rayleigh-Benard convection using the thermochromic-liquid- 
crystal technique. The number statistics and geometric properties of sheetlike/rodlike 
thermal plumes were investigated. Major findings are summarized as follows: 

(i) When observed from above, the previous-called sheetlike thermal plumes near 
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the top plate seem to be only one- dimensional structures and should be called rodlike 
plumes hereafter. These plumes evolve morphologically, i.e. convolute or spiral, to 





Figure 16. (a) Three-dimensional sheetlike/rodlike thermal plumes obtained near the 
top plate by tomographic reconstruction technique (Ra = 2.0 x 10 9 , T = 1). (b) An 
enlarged region of (a) , showing the morphological transformation of rodlike plumes into 
mushroomlike ones through convolution/spiraling. (c) An individual rodlike plume. 
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mushroomlike plumes. The width (£) of the region near the plates within which these 
rodlike plumes exist and evolve morphologically is only associated with the thermal 
boundary thickness (i.e. £ ~ several 8th)- 

(ii) The numbers of sheet like/rodlike and mushroomlike plumes, N^ eet and N£f ush , 
are found to both scale as Ra 03 . This finding suggests that the total amount of heat 
flow is dominated only by the number of thermal (rodlike or mushroomlike) plumes 
and the normalized mean heat content carried by each type of plumes is approximately 
independent on Ra. 

(iii) As the turbulent intensity is increased, large-sized sheet like/rodlike plumes are 
more easily to be fragmented into smaller ones and hence their aspect ratios decrease. 

(iv) Although sheetlike/rodlike plumes do not look like fractal objects, , power-law 
relations can be used to characterize the relationship between their perimeter and area 
and between their size and shape complexity. 

(v) Geometric measures of thermal plumes are found to have a universal 
distribution, i.e. the log-normal distribution. Thus, log-normal statistics may be used 
to model thermal plumes. In contrast, for the mixing of passive scalars measured in the 
same system, it is found that the same quantities obey log-Poisson statistics [3D] . 

(vi) The Ra-dependency of the plume area density (f r ~ i?a a23 )suggests that the 
plume number increase with Ra (N ~ Ra 3 ) can be largely attributed to the increased 
emission of plumes (as Ra is increased), rather than as a result of fragmentation of 
plumes. 
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